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We estimate exclusive density and asymmetry parameter dependent dispersion relations of vari- 
ous charged states of pions in asymmetric nuclear matter. The possibility of matter induced mixing 
of 7T° with rj is clearly exposed with the further mass modification of 7r° meson due to mixing. 
Asymmetry driven mass splitting and mixing amplitude are of the same order as the corresponding 
values in vacuum. Closed form analytic results for the mass shifts and dispersion relations with and 
without mixing are presented. Furthermore, we discuss the sensitivity of our results on the scalar 
' mean field within the framework of Quantum Hadrodynamics. 
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I. INTRODUCTION 



It is well known that the particle dispersion changes in matter due to scatterings with the medium constituents. 
This is characterized by the density(pb) dependent self-energy (£(go, |q|)) of the propagating particle. Clearly it 
depends on the energy (qo) and three momentum (|q|) of the particle. At low density T^(qo, |q|) can be calculated by 
qq ' multiplying the forward scattering amplitude with the density, which, however fails at higher density where multiple 
I/"") , scattering becomes important [1]. To incorporate these higher order effects one calculates full self-energy by evaluating 
loops at various orders. The real and imaginary parts of E(<7ch |q|) determine the in-medium mass and the decay width 
of the particle. The medium modified dispersion relations characterizing the collective excitations are obtained from 
■ the zeros of the inverse propagator involving the density dependent self-energies. 

In the present paper we plan to investigate the pion propagation in asymmetric nuclear matter (ANM). Such 
investigations are important to understand the pion(7r)-nucleon(A^) dynamics at finite density. Medium modifies the 
pion masses. Such mass shifts in nuclear matter can be used to calculate the pion-nucleus optical potential [2,3] which 
1 1 ' are different for different charged states. Furthermore, in-medium pion dispersion relations also determine the low 
mass dilepton yields in relativistic heavy ion collisions which is a subject of contemporary debate [4,5]. With this 
motivation, we here focus on the in-medium properties of pions in neutron rich matter. This is in sharp contrast with 
most of the previous calculations that deal with symmetric nuclear matter (SNM) [6,7]. The possibility of asymmetry 
driven n-i] mixing, as explored here, has not been addressed before. Such mixing also modifies the pion mass in ANM. 
It should be noted that such mixing of different isospin states and the mass splitting of various charged states in ANM 
is a generic effect and therefore interesting in itself. Furthermore, in the present study we go beyond Fermi gas model 
and incorporate N-N interactions within the framework of Quantum Hadrodynamics (QHD) [8]. 

Pion propagation in ANM is qualitatively different from SNM. This, as we show, is due to the possible splitting 
and mixing of the collective modes. Both the splitting of isospin multiplets and the coupling of the isoscalar and 
isovector oscillations are related to the ground state induced isospin symmetry breaking [9-12]. Clearly charged and 
neutral pion masses will be non-degenerate once the density dependent dressings of the pion propagators in ANM 
are considered. Such splitting may also be observed in vacuum due to virtual nuclcon-antinuclcon excitations if the 
mass(M) splitting of neutron(n) and proton(p) are taken into account (or at the quark level u-d mass difference). In 
addition under this situation tt° (pure isospin eigen state) may contain admixture of r\ meson. Such 77-77 meson mixing 
in vacuum due to the n-p mass difference has been considered in ref. [13]. But the mechanism we propose here for 
such splitting and mixing is generically different. Here, as we shall see, it is driven by the difference of the proton and 
neutron Fermi momentum (kp n ), i.e. k^ k^. This possibility, in the context of p-uj mixing, was first suggested 
by one of the present authors [11]. Moreover, to focus exclusively on the density dependent effect we neglect explicit 
symmetry breaking and take M p = M n . 

II. FORMALISM 



To describe pion-nucleon interaction we consider pseudovector coupling 
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Here = and §5 = 12.6. In Eq.l v&jv and are the nucleon and pion fields respectively and f is the isospin 
operator. The rj-NN interaction Lagrangian for pseudovector coupling can be found from Eq.l with the following 

substitution : ^- A. an d f • $ w -> Here ^- = ^ and ^- = 5.5. 
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Another essential ingredient to calculate in-medium self-energy of the tt/tj meson is the in- medium nucleon Green's 
function. In nuclear matter usual vacuum is replaced by the occupied Fermi sea which forbids on mass-shell nucleon 
propagation below the Fermi momentum because of the Pauli blocking. The relativistic nucleon propagator in this 
case is given by G(k) — G F (k) + G D (k), where the superscript F and D denotes the free and dense part respectively. 
Explicitly they are given by 

G W= k 2 -M* 2 +iC ( ) = 7* 6 (ko-e)0{k -|k|) (2) 

Here k F is the Fermi momentum and M* denotes medium modified nucleon mass. The energy of nucleon is e* = 
y/Al* 2 + k 2 . We, from now onwards, use k p and k n to denote the proton and neutron Fermi momentum respectivel. 
In QHD the effective nuclear mass is determined from the following self-consistent condition [8]. 

M * =M _ (3) 



where pf (i — p, n) represent scalar densities given by 
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The general formula for calculating self-energy is given below. 



/d 4 k 
7^34 Tr [{ iT(q)} iG p{n) (k + q){ iT(-q)} iG p[n) {k)} (5) 

Where G p („) denotes the proton (neutron) propagator and T(q) = — ^-75^ is the vertex factor. With Eq.2 and 
Eq.5 the pion self-energy due to scattering from Fermi sphere takes the following form : 

d 4 k ( f^ 2 



^^^iww T (6) 

Here T is the trace factor. For ir° 

T = Tr[ lb iG F p (k + q) l5 ftG°(k) + 75^ (k + q) 75 ^G F (k)} + [p -> n] (7) 
But for tt + (~>, fx is to be replaced by \f2fx (\/2 is the isospin factor) and 

T = Tr[ lb iG F p{n) {k + <?)75^ (p) (fc) + 75^ (n) (fc + qh^G F {p) (k)] (8) 



o P o n ^ n 




P ( \ n 

(a) 



+ n + V 

7T /^N 7T 7T 



o 



p n 
(6) (c) 

FIG. 1. (a) represents the one-loop self-energy diagram for tv°, (b) and (c) represents the same for n + and n~ respectively. 
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Now the self-eneries for tt° and n ± can be determined by evaluating the Fig. 1 (a) , 1 (b) and 1 (c) respectively. Using 
Eqs.7 and 8 in Eq.6 we get, 



S^(g) = -8 
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and 
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Where 
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with $ Pi n — 9(kp, n — |k|). The superscripts on S 0,± represent self-energies for various charged states of pion. 
We restrict ourselves in the long wave length limit i.e. when the pion momentum (|q|) is small compared to the 
Fermi momentum(fc^ n ) of the system where the many body effects manifest strongly. In this regime the concept of 
individual scattering breaks down and particle propagation can be understood in terms of the low momentum collective 
excitations of the system. In the short wave-length limit, viz. at distance scale much small compared to the inter 
particle separation (l/k F n ), one observes that the particle dispersion approaches to that of the free propagation. 
The recognition of the special role played by the distant collisions in determining the collective excitations permits 
analytical solutions of the dispersion relations [14,15]. 

III. RESULTS 

Now we present main results of our calculations. First we evaluate the pion self-energy restricting ourselves to the 
long wave length limit i.e. we consider excitations near the Fermi surface due to scattering. In effect, we assume that 
the loop momentum k ~ k F and q << k F . This permits us to neglect q A compared to the term 4(k.q) 2 from the 
denominator of both A and B in Eqs.ll and 12. Explicitly, after performing the integration we get, 
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where c = <Zo/M , v p,n — ^p,n/^ n an d e*„ = \Jk 2 n + M* 2 . In the limit of |q| = 0, Eq.13 and Eq.14 reduces to the 
following form 



S° w (go,0, kp, n ) =(ll w ql, E^(g , 0, fc p ,„) = Q, 2 nn ql ±m n 6Q,l w q , 



(16) 



where 
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^=(^—) ±\ I - i Ijj (17) 

(18) 

It is to be noted that both Q 2 .^ and tffi 2 -^ are dimensionless quantities. From Eq. 16 and 18, it is evident that the 
asymmetry driven mass splitting is an order 0(k F5 /M* 5 ) effect therefore does not appear in calculations at leading 
order in density. 

A. Isospin mode splitting of pion 



Once determined, the self-energies can be used to solve the Dyson-Schwinger equation to derive the dispersion 
relations, viz. 

ql - q 2 m 2 „ ± - S°<±( 9o , |q|) = (19) 

In ANM S° 7^ £+ ^ £~, hence the pion effective masses will split in matter. The origin of mass splitting can be 
understood by considering 7T 0,± scatterings from the Fermi surface in different channels. This would correspond to 
various cuts of Fig. 1 due to the G D term in Eq.2. It is clear that for tt^ scattering we have contributions both from 
the direct (s) and exchange (u) channels. However for the 7T + in the s channel neutron (n) density contributes while 
the u channel involves proton (p) Fermi sphere. For the ir~ self-energy the role of p and n gets reversed. While for 
7T° self-energy, both for the s and u channel, the sum of p and n densities appear. 

In the static limit (|q| = 0) following expressions for the pion effective masses (go) are found. 
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FIG. 2. The dashed, dotted and dashed-dotted lines represent the effective masses of 7r + ,7r°,7r~ respectively. Upper panel 
represents effective masses in the free Fermi gas approximation and lower panel represents the same with interaction via the 
scalar mean field. 
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From the above equation we find that the splitting is caused by the term Sd 2 .^ which for SNM (k p = k n ) vanishes 
and tt^ masses, as expected, become degenerate. 

Furthermore, it should be noted that this splitting is an higher order effect ~ k pn /e pn and therefore at the linear 
density level {k pn /e pn ) all the pion masses will remain degenerate, if the free pion masses are assumed to be so. We 
have seen that all higher order density corrections to the self-energy vanish in the static limit |q| =0. Numerical 
estimation for the effective pion masses are shown in Fig. 2 as a function of the asymmetry parameter a = (p n ~ Pp)/ Pb 
at normal nuclear matter density. Evidently in neutron rich matter the ir~ meson shows more repulsive behaviour 
than 7T° and 7r + meson which follows from Eq. 18. Moreover the 7r° mass is rather insensitive to a. We find that 
all the masses receive an overall density dependent positive corrections consistent with Eq. 17 and 18. In the upper 
panel we show the results corresponding to the free Fermi gas model, while in the lower panel effect of the interacting 
ground state is taken via Relativistic Hartree Approximation (RHA) [8]. 

The mean field, as we see from the lower panel of Fig. 2, enhances the mass splitting. It is observed that although 
the mass splitting contribute at 0(k pn /e pn ), the correction is large compared to the vacuum mass difference (Fig. 2). 
At normal nuclear matter density Am n + = -1.15(3.39) MeV, Am^o = 4.85(5.01) MeV and Am r = 11.19(6.98) 
MeV at a = 0.7(0.2) with mean field (M* = M* = 685.0 MeV). 

The dispersion relation for 7r 0± is 

(q 2 U ± ~ m* 2 ± + 7 _q 2 + (^f + ^ jj- (22) 
Effective pion masses m* ,± are given by Eq. 20. Here 

'Jinr 1 X7T7C firm 7 — ~ Ty^ (^3) 




(25) 



B. Matter induced n-rj mixing 



Another consequence of n <-» p asymmetric ground state, as mentioned before, is the isospin mode mixing. This is 
an exclusive density dependent effect. To see this we recall that 



(26) 
(27) 



The 7T° couples to p and n with opposite sign (Eq.26) while r\ couples with the same sign (Eq.27). This brings in a 
relative sign between the proton and neutron loop as shown in Fig 4 which forces mixing amplitude to vanish in SNM 
(k p = k n ) and become non-zero in ANM (k p ^ k n ). The ir°-rj mixing amplitude in the static limit is evaluated to be 



£° (<7o,0) = fi 2 q 2 



where 
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Note the difference of sign in Eq.29 and Eq.17. Clearly Q 2 n is non-zero in ANM and vanishes for e* = e* 



5 



7T 



V 



V 



7T 



n 



p 




n 



FIG. 3. This is the one-loop diagram for ir°-r] mixing 



In presence of mixing the pion and eta propagation gets coupled with each other and can be represented by 2 x 2 
matrix [8,11] 
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Shifted masses are obtained by solving the following equation. 

(q 2 - m 2 - £^( go ,0))(<z 2 - to 2 - E^(<fo,0)) - (E%(g , 0)) 2 
The Eq.31 can be written as 

4/1 a2\ 2/ *2 , *2\ , *2 *2 n 

?o(! - A ) - ?o("V +m v )+ m,„m, = 

Where 

772*0 cLud 772* are the effective masses of 7r^ and tj . Here 

ml „ nL fi* 
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Solving Eq.32 we get the modified effective masses of 7r° and rj with mixing 
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Note that for SNM, A 2 = yielding medium modified masses as given in Eq.20. 

Finally we present pionic dispersion relations in asymmetric nuclear matter with the possible ir°-r) mixing 
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FIG. 4. The dashed, dotted and dashed-dotted curves represents the dispersion characteristics for 7r + ,7r ,7r~ respectively at 
the asymmetry parameter a — 0.7. 



In SNM A 2 and S 2 vanish and Eq.36 gives the same dispersion for ir° as Eq.22. In Fig. 4 dispersion relations for 
various charged states are depicted. The pure density dependent mixing at the r\ pole is estimated to be S^ att = 
-1217.475 MeV 2 , £™ att = -1661.11 MeV 2 , at a = 0.2 and a = 0.3 respectively at normal nuclear matter density 
with coupling parameters same as that of [13] which shows that even at normal nuclear matter density the mixing 
amplitudes are of the same order as that of the vacuum mixing = —4200 MeV 2 amplitude [13]. At higher density 
the matter induced mixing significantly increases. 



IV. CONCLUSION 



To conclude, in the present work we have shown the characteristic changes of pion dispersion relations in ANM 
with the possibility of splitting and mixing of different isospin states. A density expansion of the pion self-energy due 
to scattering from Fermi sphere is made to obtain analytical results. We also have considered the modification of n° 
mass due to its mixing with the r\ meson. Furthermore, it is found that the tt~ in neutron rich matter experiences 
more repulsion than 7r 0,+ in agreement with the chiral perturbation theory calculation [3]. Both the mass splitting 
(although shown to be an 0(k^ p /e^ p ) effect) and the n-r] mixing amplitude for neutron rich system are comparable to 
their vacuum counterparts. It is found that the density dependent mixing is quite important an effect even at normal 
nuclear density for asymmetry in the range of a = 0.2 — 0.3, at higher density (and/or higher asymmetry ) both 
the mass splitting and mixing effects are dramatically higher. To focus exclusively on the density dependent effects 
the 7T-mass splitting due to n-p mass difference has purposely been neglected. Present calculation can be extended 
to include such effects which would modify the quantitative results only. We also have included scalar mean field to 
understand the effect of interacting ground state. 
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